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A novel comparison between the data and the theory is proposed for the nonmesonic 
(NM) weak decay of hypernuclei. Instead of confronting the primary decay rates, as is 
usually done, we focus attention on the effective decay rates that are straightforwardly 
related with the number of emitted particles. Proton kinetic energy spectra of iHe, \U, 
^Be, j^-B, \ 2 C, \ 3 C, and \ e O, measured by FINUDA, are evaluated theoretically. 
The Independent Particle Shell Model (IPSM) is used as the nuclear structure frame¬ 
work, while the dynamics is described by the One-Meson-Exchange (OME) potential. 
Only for the iHe, ^Li, and ^ 2 C hypernuclei is it possible to make a comparison with 
the data, since for the rest there is no published experimental information on number 
of produced hypernuclei. Considering solely the one-nucleon-induced (1AT-NM) decay 
channel, the theory reproduces correctly the shapes of all three spectra at medium and 
high energies ( E p ^ 40 MeV). Yet, it greatly overestimates their magnitudes, as well as 
the corresponding transition rates when the full OME (7r-\-K-\-r]-\-p-\-uj-\- K*) model 
is used. The agreement is much improved when only the 7r + K mesons with soft dipole 
cutoff parameters participate in the decay process. We find that the IPSM is a fair first 
order approximation to disentangle the dynamics of the 1JV-NM decay, the knowledge of 
which is indispensable to inquire about the baryon-baryon strangeness-flipping interac¬ 
tion. It is shown that the IPSM provides very useful insights regarding the determination 
the 2AT-NM decay rate. In a new analysis of the FINUDA data, we derive two results 
for this quantity with one of them close to that obtained previously. 

Keywords: Hypernuclei; Hyperon-nucleon interaction; Nuclear structure models and 
methods. 

PACS numbers: 21.80.+a, 13.75.Ev, 21.60.-n 
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1. Introduction 

The weak decay rate of a A hypernucleus can be expressed a:P 


IV = I’m + Tjvm, 


(1) 


where Tm is the decay rate for the mesonic (M) decay A —► nN, and Tnm is the 
rate for the nonmesonic (NM) decay, which can be induced either by one bound 
nucleon (1IV), T°(AN —► nN), or by two bound nucleons (2TV), T^ANN —> nNN), 
or even more bound nucleons i.e., 


r NM = r; + + ■ ■ 
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With the symbol • ■ • we indicate that additional processes, such as those induced by 
three nucleons, can contribute also. We use the superscript 0 to distinguish between 
the primary (bare) decay rates, and the effective decay rates 
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( 3 ) 


which are affected by final state interactions (FSIs) and are directly related to the 
numbers of measured single-nucleons Njv, and two-particle coincidences N„jv, as 

r^v = jv, T nN = ^N nJV , (4) 

IN w IN w 

where the number of produced hypernuclei Nw and the corresponding decay rate 
are experimentally measured quantities. The primary nucleons, in propagating 
within the nuclear environment, interact with the surrounding nucleons representing 
a complicated many-body problem generically designated as FSIs. In Eq. (J3| we 
have considered only the dominant primary decays that are later perturbed by 
FSIs, and this is the meaning of + ■ • •. 

The schematic representation of the two decay channels, when the pertinent 
dynamics is described by one-meson-exchange (OME) potentials, is shown in Fig. 
[T] This is the most frequently used model for handling the NM-decay, including 
usually the exchanges of nonstrange-mesons 7 r,p, w, and 77 , and strange-mesons K, 
and K*. It is based on the original idea of Yukawa that the NN interaction at long 
distance is due to the one-pion-exchange (OPE), with the dominant role played by 
the exchange of pion and kaon mesons. 

The OPE potential was verified quantitatively by the Nijmegen partial wave 
analysis of NN scattering in the elastic regionpl i.e., at distances larger than the 
minimal de Broglie wavelength ~ 0.5 frn corresponding to the pion pro¬ 

duction threshold. The verification of other meson exchanges is less straightforward, 
and the uncertainties in the baryon-baryon-meson (BBM) coupling constants could 
be sizeable since they are not constrained by experiments. To derive them in the 
strong sector (S vertices in Fig. [l|, the SU(3)f (flavor) symmetry is utilized. In 
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Fig. 1. Schematic representation of a) one-nucleon, and b) two-nucleon induced decays in A- 
hypernuclei when described by the interplay of weak ( W ) and strong (S') interactions through the 
exchanges of nonstrange-mesons n,p,uj, and 77 , and strange-mesons K , and K*. The W and S 
vertices are exchanged for strange mesons. 


the weak sector (W vertex in Fig. [lj, the BBM parity-violating couplings are ob¬ 
tained from the SU(6)w (weak) symmetry, while the parity conserving couplings 
are derived from a pole-model with only baryon pole resonances 00 

The BBM vertex functions also involve uncertainties in the dipole cutoff param¬ 
eters Am which, being off-shell quantities, can not be determined experimentally. 
We only know that, to have a physical meaning, they have to be of hadronic scale 
(~ 1 GeV). For instance, in different calculations, the for the off-shell pion varies 
from 0.7 to 1.3 GeV.®® In particular, the pion cutoffs of 1.2 GeV, and 0.8 GeV 
were used to describe the electromagnetically induced two-nucleon emission pro¬ 
cesses ( 7 ,pn) and ( 7 ,pp)p which are quite similar to the 1AT-NM decays (A p,np), 
and (A n,nn). We use here the Am’s from Ref.,® and also those proposed in Ref00 
to account for the transition rates T p and r„ in the s-shell hypernuclei. Q 

The short-range correlations (SRCs) between the emitted nucleons nN, and 
nNN can also affect significantly the BBM couplings. Parreno and Ramos have 
shown that they can diminish the value of T ) 1 by more than a factor of two . 11 
Nothing has been stated so far regarding the effect of the SRCs on the 2JV-NM decay. 
The theoretical scene becomes even more complex when effects of the quark degrees 
of freedom , 12 13 the 27r-exchangesp®£ 17 and the axial-vector ai-meson exchangdUE 17 
are considered. 

The M and 17V-NM decays have been observed experimentally in the pioneering 


a The dependence of the NMWD transition rates on the values of Am, and therefore on the BBM 
vertex functions, is thoroughly discussed in Table 1 and Fig. 1 of Ref® 
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measurement performed more than 50 years ago by Ruderman and R. KarplusP^ 
Conversely, the experimental observation of the 27V-NM decay, which was predicted 
by Alberico et al , 19 in 1991 (see also Ref! 20 }), has been reported only in recent 
years at KEKpH and at FINUDA. 22 ^ Both groups announced a branching ratio 
r°/r nm ~ 25 — 30 %. The first group obtained this result from the single and 
double coincidence nucleon spectra in ^ 2 C, and the second one from proton kinetic 
energy spectra in ((He, \ Li, 9 Be, ” B, ^ 2 C, \ 3 C, ^ 5 N and ^ 6 0. A branching ratio for 
the 21V-NM decay channel of this magnitude is consistent with the prediction made 
by Bauer and GarbarinoP^ On the other hand, the one-particle proton and neutron 
kinetic energy spectra in ((He measured at BNlP^are accounted for reasonably well 
theoretically by considering only the 17V-NM decay modelSEI 

The above mentioned experiments, together with several others performed dur¬ 
ing the last few decades(HS 27 MHrepresent very important advances in our knowledge 
about NM decay. Explicitly, these advances are: 1) new high quality measurements 
of the number of single-nucleons Njy, as a function of the one-nucleon energy En, 
and 2) first measurements of the number of two-particle coincidences N n jv, as a 
function of: i) the sum of the kinetic energies E n jv = E n + En, h) the opening 
angle 9 u n , and iii) the center of mass (c.m.) momentum P h n = |p„ + p N |. On 
the theoretical side this implies a new challenge for nuclear models which have 
to explain, not only the IN- and 2iV-NM decay rates, but also the shapes and 
magnitudes of all these spectra, testing in this way both the kinematics and the 
dynamics. 

Recently, Bauer, Garbarino, Parreno and RamoP® have obtained good agree¬ 
ment with KEK data}—] considering both the one- and the two-nucleon induced 
decays in the framework of the Fermi Gas Model (FGM). These authors have also 
analyzed the proton kinetic energy spectrum in \ 2 C measured at FINUDA J^but 
no theoretical analysis of the remaining spectra has been done so far. In the present 
work, we present for the first time the calculation of all proton kinetic energy spectra 
studied in the above mentioned experiment. 

Since i) the NMWD is dominated by the 1IV-NM decay, and ii) the 2iV-NM 
processes and the FSIs contribute mainly at low energy, it is reasonable and useful 
to compare the experimental spectra with theoretical calculation when only the 1 N- 
NM decay mode is considered.^] It is obvious that there will be no agreement at low 
energies between the experimental and theoretical spectra. But this disagreement 
is not crucially important, since we are basically interested in disentangling the 
strangeness-flipping interaction among baryonsPOf course, both the 2iV-NM decay 
and the FSIs, as well as the SRCs, are interesting physical phenomena in themselves, 


b Such a comparison is analogous to those done between the experimental data on electron- 
nucleus and charged-current neutrino-nucleus scatterings that include the FSIs, with the plane- 
wave impulse approximation which doesn’t include these processes I2£|E2] Moreover, the results of 
our analysis are fully robust, in the sense that they will be valid even after the inclusion of the 
FSIs and the 2/V-NM decay. 
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but they teach us little about the basic nonmesonic decay. Moreover, as indicated in 
([3]) they can not be treated separately, and it is not known whether they contribute 
coherently or incoherently. That is, it can even happen that they partially cancel out 
(for instance, T^ FSI and terms in ©), as do the divergences and the vertex 

corrections in the QED, because of Ward identity. More specifically, and as already 
pointed out in Ref.,® the issue of FSIs in the NMWD is a tough nut to crack, 
and there is no theoretical work in the literature encompassing all aspects of these 
processes. For this reason, before having a reasonable control over all the physics 
that they involve, it may be useful, or even preferable, to discuss the experimental 
data without the FSIs. This is what we do here. 

The content of this article is as follows. Our method to compare the experimental 
data with theory for the NMWD is explained in detail in Sec. [2] The main formulas 
used to calculate the proton spectra for the liV-NM decay within the Independent- 
Particle Shell Model (IPSM) is presented in Sec. [3] The parameterizations that 
are used for vertices are listed in Sec. [3] also. The only novelty here is that the 
proposed BBM vertex functions are rarely used in the literature. The calculated 
spectra for all hypernuclei are presented in Sec. [4j where we make a comparison 
between theory and the FINUDA data for ((He, ((Li, and ^ 2 C. The extraction of r° 
from experimental spectra is reanalyzed in Sec.[5j Finally, Sec. [G] presents the final 
remarks and conclusions. 


2. Relationship between Experiment and Theory 

While Fjv and r n jv are experimentally observable quantities, the bare decay rates 
F°, r°, and r° are not and have to be derived from the data, employing different 
extraction procedures which frequently involve approximations that are question¬ 
able. (One example will be illustrated here.) Moreover, the primary decay rates are 
ill defined, since they depend on the model that is used to describe the nuclear 
structure of the hypernucleus, as well as on the SRC, etc. 

Note that we include the FSIs in the definition of r^r, and r ra j\r which is not 
commonly done in the study of the NMWD EM?' But, there are other processes 
in nuclear physics where the FSIs participate in the definition of the decay rates. 
The best known phenomenon is nuclear /3-decay, where the transition rate depends 
on the FSIs caused by the Coulomb attraction of the emitted electrons. (See, for 
instance, Eq. (5.11) in Ref. 41 where the FSIs effects are approximated by the Fermi 
function.) The main difference between the FSIs in the leptonic and nonmesonic 
weak-decays is that while in the first case they are easily evaluated, in the second 
case they are very complicated-® and beyond the scope of this work. 

The FSIs are usually simulated by a semi-classical model, developed by Ramos et 
aZ.j^ancl denominated Intranuclear Cascade (INC) code. This code interrelates the 
primary rates with measured rates ©>• More recently, the FSIs were evaluated 
with a time-dependent multicollisional Monte Carlo cascade schemeJ2?EI imple¬ 
mented within the CRISP code (Collaboration Rio-Sao Paulo), which describes, 
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in a phenomenological way, both the nucleon-nucleon scattering inside the nucleus 
and the escape of nucleons from the nuclear surfaceP^MHl The CRISP code, as all 
INC codes, is tailored to simulate the experimental data, and not for describing 
theoretically the FSIs. As such, it involves a normalization procedure for the ex¬ 
perimental data, which washes out all of the information on the NMWD dynamics. 
More specifically, the spectra of rj evaluated from a Shell Model (SM) are the main 
ingredients for establishing the initial condition to start the CRISP cascade process 
and to calculate in this way the FSIs. (The primary spectra of r° should also be in¬ 
cluded within the initial conditions in a more complete model, but we do not know 
yet how are they evaluated within the SM.) Yet, because of the normalization, the 
spectrum perturbed by the FSIs turns out to be the same for different primary 
spectra. (One obtains the same V p ’ FSI spectra for different spectra, etc.) We 
were not able, so far, to find out how to circumvent this problem of normalization. 
Moreover, not all FSIs are considered within the INC codes. Which additional FSIs 
contribute to the NMWD spectra and decay rates, and how and which of them 
should be included in the calculation are nontrivial questions. Some candidates are 
discussed in RefP^ 

The information on the dynamics also is lost when the decay rates Tjv, and 
r n N are normalized to T^m- This normalization is the usual procedure; see for 
instance 36 Eq. (7)], and |49 Eqs. (14),(15)]. With this normalization, the spectra 


depend on the phase space and the FSIs, but very weakly on the NMWD dynamics. 
The same happens when the transition density is normalized to the decay rate 
(See [57| Fig. 3].) In contrast, as explained below, we do not normalize any of 
the calculated transition rates, and this allows us to inquire more deeply into the 
NMWD mechanism. [] 

The measurement implies the counting of the numbers of emitted protons AN£, 
and the errors <5ANp, corrected by the detection efficiency, within m energy bins of 
width A E. The total number of emitted protons and the resulting errors are 


N P = ANp, dN p = 

<=1 \ 


5> AN p) 2 ’ 


( 5 ) 


i— 1 


where the summation goes over proton energies E l p larger than a given threshold 
energy E thres . 


c It might be useful to draw a parallel with electromagnetic decay. When the electric E2 and the 
magnetic Ml multipoles are the lowest allowed transitions, both may contribute significantly to 
the total rate T 7 = Y 1 {E2) +r 7 (Ml), since the electric transition may be enhanced substantially 
above the single particle estimate due to collective effects. The comparison between theory and 
data is always done for the decay rates T 1 (E2) and T 7 (M 1), separately. There is no physical 
motivation for comparing the ratios r 7 (.E2)/r 7 and r 7 (Ml)/T 7 , since they are less sensitive to 
the nuclear structure effects than the individual decay rates. 
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The corresponding total decay rate with its error read 


r p = E AI l>’ = 


i— 1 




\ i=l 


( 6 ) 


where, as seen from 0, the decay rates AT l (E p ) with errors JAr* are given by 


Ar i(E p ) = ^AN ;(E P ), 


and 


<5Ar; = *Ean* 
p Nw 


/ (5AN* 
V AN* 


( ST W 
V 


2 (5N W ' 2 

\N W 


1/2 


( 7 ) 


( 8 ) 


with and 5N\y being, respectively, the experimental errors on r^, and Nw- 

Thus, to evaluate the experimental decay rates we need to know the values of Tw 
and Nw for each hypernucleus. As usually done all T’s will be given in units of Ta, 
the total decay width of the free A. 

For the first quantity we can use the relationship 


T W {A) = (0.990 ± 0.094) + (0.018 ± 0.010) A, (9) 

which was derived in RefP^ from a linear fit to the known values of all measured 
hypernuclei in the mass range A = 4—12. 

But unfortunately there is no experimental information about Nw in the liter¬ 
ature. Only the ratio 


R p = 


N p 

Nw’ 


( 10 ) 


for the eight hypernuclei discussed here were presented at a conference PAl for the 
threshold energy E thres = 15 MeV. We will use, however, only the results for )\He, 
^Li, and \ 2 C, since only these were published in a refereed physics journal so farP2l 
We list them in Table [lj together with values of N p from Ref. and the resulting 
estimates for Nw from (101. The proton decay rates, evaluated from 


T p = R p Tw > (11) 

are also shown. Since, according to ([9]), the value of Fw is close to unity, the last 
result implies that the ratio R p is basically the proton decay rate. 

The theoretical analogs of |6]) and ([7| are, respectively 

Tl h = J S p (E p )dE p , (12) 

and 

AF l h (E p ) = S P (E P )AE, (13) 

where the spectral function S p (E p ) depends on the theory that is used to evaluate 
the NMWD, which not yet has been discussed. 
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Table 1. Values of R p from Ref!^ and N p from Ref!^ for Eth res = 15 MeV, together with the 
resulting estimates for Nw, and F p = R p Tw- 


Hypernucleus 

R p 

N p 

Ntv 

r P 

iHe 

0.25 ±0.07 

262 ± 25 

1047 ± 391 

0.27 ±0.11 

iLi 

0.37 ±0.09 

259 ± 21 

700 ± 226 

0.41 ±0.16 

12p 

A ^ 

0.43 ±0.07 

678 ± 38 

1576 ± 344 

0.52 ±0.18 


Instead of comparing the experimental transition rates with the calculated rates, 
we can compare directly the number of measured protons with the calculated quan- 
tity 

A<0E P ) = ^S(E P )AE, (14) 

1 w 

where Nw /Fw is just a proportionality factor. 

All the above is completely general. The theory that is used can be as compli¬ 
cated as necessary to properly interpret experimental data. But it can also be very 
simple and still lead us to correct conclusions about the underlying physics. |^] Such 
a model is described below. 

3. Independent Particle Shell Model for the Spectral Function 

The IPSM has been used for more than twenty years in the evaluation of the 1N- 
NM decay ratesbut only in recent years was it applied for the description of 
different spectral densities Sn{E), S h n(E), (cos 0) , and 5„jv(P) 0OUEIllAUEUE23 
We briefly sketch here the main assumptions that are made in this model, and give 
the resulting theoretical expression for the proton kinetic energy spectrum. 

The assumptions are: (i) the initial hypernuclear state is taken as a hyperon A 
in a single-particle state j\ = 1 Si/ 2 weakly coupled to an (A — 1) nuclear core of 
spin Jc, i.e., | Jj) = \(JcjA)Ji)', (h) the nucleon (N = p,n ) inducing the decay is 
in the single-particle state j n {j = nlj); (iii) the final residual nuclear states are: 

| Jf) = \(Jc3n 1 )Jf)\ (iv) the liberated energy is 

= A + e\ + Sjy, (15) 

where A = Ma — M p = 177.33 MeV, and the e’s are experimental single-particle 


d The model used to describe a process should, in principle, describe all the involved physics. This 
is a desirable condition, but it is not in any way essential. Useful models are those that allow us to 
infer consequences consistent with the observations. More precisely, a model is a simplified version 
of the process, and the model designer decides which features to consider. 
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energies (s.p.e.) Q , and (v) the c.m. momenta and relative momenta of the emitted 
particles are: 


Pun = V ( A - 2 )(2M,- Pi ~ P 2 n ), 


Pn N -J MA N 4(^4 _ 2) P ™ N ' 


It follows that the liV-NM decay rate is given bj m 

= ^ ^ 1 "'/y ? = I P"]y(PnN i PnN^dQnN : 

j 

where dfl nA r is the phase space factor, and 

J=j+ 1/2 


P N (P,P) = 


1 + 5 


Nn 


E F NJ E \M{plPLXSJT-j A ] N JT)\ 2 , 


J=\j-l/2\ SILXT 


with 


M (plPLXSJT ; j\j N JT) = E [i _ (-)*+ s + T ] 


x 0 L (P)(lLXSJT\V(p)\j A j N JT), 


and 

0l(P) = J R 2 dRj L (PR)R 0L (b/V2,R). 

Here L 1 l , and A are, respectively, the c.m., relative, and total orbital angular 
momenta (A = L + I), while V is the transition potential, and b is the harmonic 
oscillator length parameter. 

The kinematics of different spectra Sn depend on dfl nA r and on the over¬ 
lap 0l{Pun), while the decay dynamics is contained in the matrix element 
(lLXSJT\V(p n N)\jAjNJT). The information on nuclear structure is enclosed in 
the spectroscopic factors F° NJ , which account for the Pauli Principle within each 
single-particle shell jpj. In the general case, they are given by 

*1/ = (2./ + i)E{ J j f N 3 j F ] IM4JI j £)I 2 , ( 16 ) 

where (Jcll a ] lV II^F) are fractional parentage coefficients (FPCs), and the sum¬ 
mation goes over the n final states Jp in the residual (A — 2) nuclei, with the same 


e The schematically drawn energies in [59[ Fig. 7] are the experimental s.p.e., which can be identified 
with the SM s.p.e. only in closed shell nuclei, such as 16 O. For open shell nuclei, which is the case 
of 12 C, the experimental s.p.e. are frequently identified with the quasiparticle energies, which 
include the effect of pairing correlations, and could be quite different from the SM s.p.e.. This was 
done, for instance, in 60] Table IV], where the experimental and SM s.p.e. are listed, respectively, 
in columns two and five. Moreover, the correct p 3 / 2 > P 1/2 and ^ 5/2 experimental energies read, 
respectively, —15.96, —1.95, and 1.61 MeV P^l 
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spin, parity, and isospin, and different excitation energies. Such a detailed descrip¬ 
tion could be redundant since, in the NMWD, we evaluate the inclusive decay rate, 
without being interested in exclusive processes that feed each of the individual final 
states Jp. Thus, within the IPSM, the spectroscopic factor becomes much simpler 
since | Jp) —> \(JcJn 1 )Jf), and the summation goes only over the values of Jf that 
fulfill the constraint ,Jc — Jn < Jf < Jc + Jn- The values for Jj and Jc are taken 
from experimental data and, for the hypernuclei of interest here, are listed in Table 
I of RefP 9 The resulting factors F 3 * NJ are listed in Table II of the same paper. 

The spectra are obtained from the differentiation of Tat with respect to Eat, 
cos 0 n Ni E n at, and P n N- In particular, for the kinetic energy spectrum, one has: 

S N (E N ) = Y J S j N {E N ), (17) 

3 


with 


where 


and 


Finally, 


o M 3 /*+l / fP 

S 3 n (E n ) = (A- 2)- / dcos9 nN J-=^E n l jN (pP) 

71 J -1 V h N 


E' n = (A- 2 )(A - l)Aj vr - E n [(A - l) 2 - cos 2 e nN ] 


E n — 


E' n - t/En cos 0 nN 


1 2 


(A-i)- 2 . 


r N = J2 r N = E [ QN S J N (E N )dE N , Q% = J A^, 


(18) 


with Q J N being the single-particle Q-values. 

The outline of the numerical calculation is the following: 


(1) The transition potential F(p n jv) for the emission of the nN pair, contained in 
d~j\[j l (PnN )i is described by three OME models, namely: PI - The full pseu¬ 
doscalar (Tr,K,rj) and vector ( p,u,K*) meson octets (PSVE), with the weak 
coupling constants, and dipole form-factor cutoffs A m from Refs. pus P2 - 
Only one-(7r + K) exchanges (PKE) are considered, with the same parametriza- 
tion as in the previous case, i.e., with cutoffs = 1.3 GeV and Ap = 1.2 GeV 
fromjSand P3 - The soft n+K exchange (SPKE) potential with cutoffs A w = 0.7 
GeV and A K = 0.9 GeV fromHE 

(2) The SRCs acting on final nN states are incorporated phenomenologically 

through Jastrow-like SRC functions, as used within both finite nuclei calcu¬ 

lations PHHH and Fermi Gas Model (FGM) calculations GHIHI 
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4. Proton Decay Rates and FINUDA data 

The calculated transition densities AT^^Ep), for A E = 10 MeV and evaluated 


from (13) with S p (E p ) given by are shown in Fig. [2] for the hypernuclei mea¬ 
sured by FINUDAP As expected, in all cases the spectra strongly depend on the 
parameterization that is used for the transition potential, while their shapes are 
the same for all practical purposes. The experimental values of AT p (E p ) for ^He, 
A Li, and \ 2 C, evaluated from (J?]) for the values of AN P (E P ) shown in [24] Fig. 1], 
are also displayed in the same figure. Their errors were calculated from (18]). 




0 20 40 60 80 100 120 140 0 20 40 60 80 100 120 140 

Kinetic Energy (MeV) Kinetic Energy (MeV) 

Fig. 2. (Color online) Calculated proton kinetic energy spectra Ar p h {E p ) for three different tran¬ 
sition potentials specified in the previous section. For ')(He, ^ l.i, and ^ 2 C are also shown the values 
of Ar p (E p ) evaluated from 0 that correspond to the AN p (E p ) shown in [24[ Fig. 1]. 


With the parametrization PI, the theory greatly overestimates the experimental 
spectra at medium and high energies ( E p > 40 MeV), underestimating them at low 
energies. The discrepancy can not be settled by simply opening a new 2A-NM 
decay channel induced by two nucleons, since this decay mode, although capable 
of producing additional particles at low energies, is unable to lower the transition 
strength at high energies. Nor is it possible for the FSIs to solve the problem, since 
they can hardly change the total transition density induced by a proton. They 
can only remove a portion of the strength from high energy to low energy. It is 
self-evident from Fig. [2] that such a mechanism cannot be successful in the present 
case. 

Improved agreement is obtained in the P2 model, which means that the incor- 
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poration of vector mesons, instead of improving the agreement, makes it poorer. 
However, the high energy part of the ](He spectrum is reproduced fairly well only 
with the parametrization P3. The decrease in magnitude of proton spectra in going 
from PI to P2 is due to the well known fact that the parity-violating contribu¬ 
tion of the vector mesons to the proton transition rates is quite sizable (see, for 
instance, 53, Table IV]). On the other hand, the strong variation of the same ob¬ 


servable in ^He with regard to the BBM vertex functions is thoroughly discussed 
in |8) Fig. 1]. 


Table 2. Transition rates F), calculated from (J 1 2|) for all three parametrizations and with two 
different threshold energies: Ethres = 15 MeV, and Ethres = 40 MeV. For )He, ^Li, and C the 
values of T p and their errors, evaluated from are also shown. 


Hypernucleus 

Bp 

p PI 

P 

pP2 

1 p 

pP3 

V 

Ethres = 15 MeV 





a He 

0.27 ±0.05 

0.466 

0.360 

0.185 

A Li 

0.41 ±0.07 

0.531 

0.415 

0.228 

iBe 


0.627 

0.494 

0.275 

a b 


0.667 

0.527 

0.289 

12 

A ^ 

0.52 ±0.07 

0.792 

0.627 

0.343 

13p 

A ^ 


0.776 

0.614 

0.336 

a 5 N 


0.821 

0.651 

0.365 

16 o 

A U 


0.906 

0.718 

0.399 

Ethres = 40 MeV 





a He 

0.19 ±0.04 

0.428 

0.331 

0.171 

A.Li 

0.27 ±0.05 

0.488 

0.381 

0.204 

9 Be 


0.555 

0.437 

0.243 

aB 


0.607 

0.481 

0.263 

12p< 

A ^ 

0.37 ±0.06 

0.719 

0.569 

0.311 

13p 

A ^ 


0.706 

0.559 

0.305 

a 5 N 


0.721 

0.573 

0.320 

a 6 0 


0.801 

0.635 

0.352 


At this stage, it might be useful to recall that, in all theoretical descriptions 
of the NMWD, only three body final states have been considered, implying that 
the residual nuclei are necessarily bound. Obviously, this is done for simplicity, but 
does not always occur. The most emblematic case is that of ^B, which has been 
considered in many theoretical studies done so far. However, its parent nucleus in 
the neutron channel, n B, is unstable and disintegrates into p+2a with a half-life of 
8 • 10~ 19 s, which is very short when compared with the half-life of B. Among the 
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nonmesonic decays analysed here, the same occurs with the 5 He nucleus, which is 
the residual nucleus for the proton NMWD of ^Li. In fact, it is unstable to particle 
emission, decaying into p + 4 He with a half life of 70(3) • 10 _25 s. This time is much 
shorter than the lifetime of \Li and, therefore, the instability of 5 He could be the 
cause of the discrepancy between theory and experiment at high energy. 


Table 3. Spectroscopic factors and the transition rates T v for )(Li evaluated in the intermediate 
coupling model with the wave function amplitudes agg given in Ref. EH which are also listed. 


a-LS 


LS 

fitffl 

fit5 

CK 616 

CK6 16 

CKPOT 

01 

0.9873 

0.9906 

0.9576 

0.9484 

0.9847 

12 

-0.0422 

-0.0437 

-0.2777 

-0.3093 

-0.1600 

10 

-0.1532 

-0.1298 

-0.0761 

-0.0703 

-0.0685 

F j 

r NJ 






JJj 






0,Pl/2 

0.0320 

0.0311 

0.0693 

0.0790 

0.0419 

l.Pl/2 

0.3023 

0.3032 

0.3026 

0.3022 

0.3042 

4j P 3/2 

0.5255 

0.5259 

0.4973 

0.4903 

0.5175 

2,P3/2 

0.1403 

0.1399 

0.1307 

0.1286 

0.1364 

Tp 

0.2444 

0.2444 

0.2475 

0.2483 

0.2453 


Similarly to what was done in Figure [2] for proton kinetic energy spectra 
AT th {E p ), the theoretical results for the proton decay rates (18) for two differ¬ 
ent threshold energies ( E t hres = 15 MeV, and E t hres = 40 MeV) are shown in 
Table [ 2 ] For |He, ^Li, and \ 2 C the experimental values of and their errors, 
evaluated from ([6]), are also shown. 

All that was previously stated when comparing the theory with the data in 
Figure [2] also applies here. In particular, since the theory does not include the 
2N-NM channel, the experimental transition rates must always be larger than the 
theoretical rates for Ethres = 15 MeV. This condition is satisfied only for the 
parameterization P3, suggesting that the other two sets of parameters would not 
be appropriate. On the other hand, as the FSIs, which remove the density transition 
from the high energy region, are omitted in the calculations, all three should 
be larger than for E t h re s = 40 MeV. From Table [2j we see that within the 
experimental errors this is indeed the case. However, in calculations P2 and P3, the 
differences between the data and the theory are too large to be entirely attributed 
to the lack of the FSIs in the theory. 

It would also be interesting to compare the experimental results shown in the 
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upper part of Table [|] with the calculation done by Itonaga and Motoba 17 with a 
more elaborate SM than that used here, employing the Tr + 2'ir/p + 2'K/a + ui + K + 
pn/ai + crn/ai exchange potential. They obtain: r°(^He)= 0.237, and T p (\ 2 C)= 
0.534. In the first case, T p is greater than T p as it should be, but it is lower in the 
second case which is not correct. 




0 20 40 60 80 100 120 140 0 20 40 60 80 100 120 140 

Kinetic Energy (MeV) Kinetic Energy (MeV) 

Fig. 3. (Color online) Experimental datJSI for proton kinetic energy spectra for the number 
of particles AN P (E P ) are compared with the IPSM results for three different OME potentials. 
Theoretical results have been evaluated from equation Also shown are the Gaussian-function 
fits of each proton spectrum from 80 MeV onwards, performed in RefPil 


It can be argued that a derivation of the F j NJi based on the jj single-particle 
model, is not fully appropriate for light nuclei, such as Li and Be, with the in¬ 
termediate coupling model preferred over pure jj coupling for the core nuclei. In 
fact, the structure of these nuclei is closer to LS coupling than to an assembly of 

Table 5] where the 


63 


P 3/2 valence nucleons, as can be seen, for instance, from 
core state Jc = 1 + in ^Li is basically a pure 3 S state. In this case, instead of the 
spectroscopic factors F^( 2 =5/8, and F^ 2 =5/8 given in Ref. 59 Table 2], one 


jvi — u i u > ± N2 

has F p n 3 ( 2 = 19/36, F^ 2 = 5/36, F^ 2 = 1/36, and F^{ 2 = 11/36. Within the 
parametrization P3, the last spectroscopic factors yield T p = 0.244 which is only 
slightly larger that the value 0.228 shown in Table[l] However, to justify even more 
reliably the jj coupling, we have recalculated T p with five different wave functions 
evaluated in the intermediate-coupling model, and cited inP^l Their amplitudes a^s 
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are listed in Table[3] The spectroscopic factors ([6| are evaluated from the expression 


= 3(2J+l)(2j p + l) 

E 

(2 J F + 1) 

Jf 

= 1/2,3/2 

x 5^(2L + 1)(2S+1) 


f 1 ih y 

ClLS 

1 1 

LS 


i JJ 


1 1 I 
2 2 

J jp Jf 


(19) 


The results are shown in Table [3j from which it can be concluded that the difference 
between the jj coupling and the intermediate coupling is < 10%. The physical 
reason for this fact is that T p is an inclusive quantity, so it is not acutely relevant 
if the proton decays from orbital p 3 / 2 or Pi/ 2 . 

We would like to stress that the way to compare the theory with data as done 
here, as well as in a previous paper , 9 Ellis conceptually different from the traditional 


way) 


This can be seen immediately by confronting expression (13) 


with 36 Eq. (7)]. Instead of comparing different bare proton contributions T 9 , T 9 , 
T°, • • •, which are not directly measured but are extracted by the experimentalist 
from the data, we compare the total decay rate T p and the corresponding spectra, 
which include all of the protons that come from the NMWD. Different extraction 
procedures are not unambiguous, as we have discussed in the first version of the 
present worlPH regarding the derivation of T 9 by both KEKpJ and FINUDA.^1 
More details on the second extraction procedure are given in the next section. 

Fig-! is just the replica of Fig. [2] except for the numbers of protons, with 
AN ^(Ep) evaluated from (14). For reasons of completeness, in the last figure we 


show the spectra of all measured hypernuclei. Obviously, the two figures lead to 
the same conclusions. The Gaussian-function fits of each proton spectrum from 80 
MeV onwards, that were performed in Ref.^ and that will be discussed in the next 
section, are displayed also in the last figure. 


5. Extraction of branching ratio r^/V° NM from the data 

The determination of the branching ratio F 9 /r ° NM by FINUDA^ is based on the 
partition of the total number of detected protons N p into low, and high energy 
regions populated, respectively, by N p = N p (E p < E part ), and N p = N p {E p > 
-EUrt) protons relative to the partition energy E part . After assuming that all 2N- 
NM protons are contained within N p , they define the ratio 

= n 9 < + n 9 + Np SI r 9 </r 9 + rg/rg + n< si /n° p 
N p N 9 + N 9 + N F si 1 + r 9 /r 9 + n fsi /no 5 1 J 

where N 9 and N 9 are, respectively, the numbers of protons induced by the 1N-NM, 
and 2N-NM decays, and N F $i = N p SI + N^ SJ is the total number of particles 
produced by the FSIs. 

The next steps done in Refill are not supported by sufficiently firm physical 
arguments. Namely, it is assumed: i) that the proton spectra from 80 MeV onwards 
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are due entirely to protons coming from the A p —t np reaction, ii) that they can 
be fit by Gaussian curves shown in Fig. [3j and iii) that the maxima of these curves 
correspond to the partition energies £ part . All of this yields 

n< _ o.5 + rj)/r° + n < si / n° p 

Np 1 + r°/r° + n FS i/ no ’ 1 > 

since the Gaussian curves are bell shaped, satisfying always the condition F^ = 
r p(E p < £ max ) = F°> = r°(Ep > E max ). The resulting partition energies E™ UDA 
are listed in the second column of Table |U 


Table 4. Mean values of FINUDA Gaussian fitd^(in units of MeV) are confronted with the maxima 
of proton spectra with the theoretical maxima (third column), as well as with the energies -Seven 
for which the proton strength is evenly distributed, i.e., (fourth column). 


Hypernucleus 

^FINUDA 


E ltn 

A He 

68.5 ±4.1 

77.5 

75.0 

A.LI 

76.7 ±5.2 

73.5 

72.0 

^Be 

78.2 ±6.2 

69.0 

69.0 

a b 

75.1 ±5.0 

69.0 

70.5 

12p 

A ^ 

80.2 ±2.1 

67.5 

70.5 

13p 

A ^ 

83.9 ± 12.8 

67.5 

70.5 

a n 

88.1 ±6.2 

61.5 

66.0 

A° 

93.1 ±6.2 

61.5 

66.0 


Next, FINUDA approximated (21) by a linear function of the mass number A, 


R(A ) = a + bA , 


( 22 ) 


where 


_ o.5 + r°/r° 
i±r°/r° 


(23) 


does not depend on A. Finally, a y 2 fit for R(A) was done for the energies F™ UDA 
to obtain the values of a and b that are shown in row A of Table [5j together with 
the resulting r^/Tp, and rD/T)^ derived from 

a- 0.5 F° a - 0.5 

rf~ l-o ’ ~ (i-a)r°/r°±o.5’ 

for the experimental value T° /F° = 0.48 ± 0.08, measured by KEKpSl 


(24) 
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The FINUDA procedure to extract the value of from a series of kinetic 
energy spectra, by separating them into low and high energy regions, looks phys¬ 
ically sound. However, we shall soon see that it is very sensitive to the separation 
procedure. On the other hand, in the fitting of the spectra with Gaussian curves, 
FINUDA implicitly assumes the absence of Np SI , which is not only unrealistic, but 
also not necessary. 

Before separating the spectra, we compare the calculated energy locations £U* ax 
of the spectra maxima with the maxima of the FINUDA Gaussian fits E™ UDA . 
From Figs. [2] and [3j one immediately notices sizeable differences. Moreover, from 
Table |4l one sees that, while the maxima £ FINUDA increase from 68.5 MeV to 93.1 

I I max 

MeV in going from ((He to (( 6 0, the energies FU* ax decrease from 77.6 MeV to 61.5 
MeV. On the other hand, while the Gaussian curves are bell shaped, the theoretical 
1N-NM spectra deviate significantly from a symmetrical shape. More precisely, the 
calculated F^ is greater than F^ in ((He, ((Li, and ((Be, and smaller for the 
remaining hypernuclei. 

The reason for this can be understood from an inspection of Fig. |4j where the 
spectra of ((He and \ 6 0 without recoil (a) and with recoil (b) are displayed. First, as 
expected, the recoil effect is sizeable in ((He where only the orbital s 1 / 2 contributes. 
Second, the ((He spectrum does not have the symmetric bell shape, mainly because 
the single kinetic energy reaches its maximum value rather abruptly at ~ 115 
MeV due to the recoil factor (A — 2 )/(A — 1) = 3/4 in Eq. (18 1 for Qp 1/2 ', this 
effect, however, does not modify the value of E max = 77.6 MeV, but causes to 
be appreciably larger than F/^ (F^ = 0.55 r°). Third, in the case of \ 6 0, three 
partial waves (s 1 / 2 , P 3 / 2 , andpi/ 2 ) contribute with different heights and widths, the 
convolution of which is a nonsymmetric proton spectrum with T p < = 0.43 here 
the energy E max becomes significantly smaller because the energy Ap 1/2 , given by 
(15), is ~ 25 MeV smaller in \ 6 0 than in ((He. Briefly, as the value of A increases, 


the average value of the binding energies £a and e 3 N turns out to be larger (see 
Ref. 66 Fig. 11]), which makes the energy position of the maximum E^ gx for the 


liV-NM proton kinetic energy spectrum become increasingly smaller. The numerical 
results for £7/ x are shown in TableEJ where they are compared with the peaks of the 
FINUDA Gaussian fits E FINUDA f2SNote that E th decreases with the mass number 

max max 

while E™ UDA increases. Thus, they differ from one another quite significantly and 
the difference increases from 9 MeV in ((He up to 32 MeV in / 6 0. 

The energies E max are closely related to the liberated energies, i.e., to the Up¬ 
values since the latter should, in principle, also decrease when the A J p decrease. 
However, for light hypernuclei, this decrease is largely offset by the recoil effect, as 
shown by Eq. (18). The final results, displayed in Figs. [ 2 ] and [3j demonstrate that 
the Q-values are roughly constant and within the energy interval of ~ 115 — 135 
MeV, which is consistent with the data within experimental errors. [(] 


1 One should also mention that it is assumed here that the residual nucleus is emitted in the 
ground state, and that consideration of excitation energies could further diminish the Q-values. 
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Fig. 4. (Color online) and ^ e O spectra: (a) without recoil, and (b) with recoil, evaluated 
within the IPSM for the P2 OME transition potential. Contributions of different orbitals S 1 / 2 , 
P 3 / 2 , and P 1/2 to the total ^ e O spectra are also displayed. 


Table 5. Results for the parameters a and b, and the corresponding ratios T^/Tp, and r^/r^M 
for FINUDA datapH and different partition energies for N<, and N>: A) derived in RefP®l with 
^max UDA , B) and C) obtained here with E^ ax , and E ^ en , respectively. 


case 

a 

b 


r °/ r NM 

A 

0.654 ±0.138 

0.009 ±0.013 

0.43 ±0.25 

0.24 ± 0.10 

B 

0.711 ±0.101 

-0.010 ±0.008 

0.73 ±0.61 

0.33 ±0.10 

C 

0.656 ±0.101 

-0.004 ±0.008 

0.45 ± 0.43 

0.23 ±0.09 


In view of the differences between the FINUDA spectra and those calculated 
here, one is immediately tempted to repeat the above analysis using the energies 
E th instead of _£; FINUDA . This was done and the results for R are shown in the left 

max max 


panel of Figure [5J along with the corresponding linear \ 2 fit. The large difference 
between the two sets of maxima gives rise to large differences between the values 
of the ratios R in Ref. [24j Fig. 2] and those derived here. The parameters a and 
b, and the ratios T^/Tp, and F!j/r°- M obtained in this way are listed in row B of 
Table[5j The value of a is not very different from the previous case but, as expected, 
b is now negative, and the resulting is significantly different due to its strong 
sensitivity on a in (24). 

One must not forget here that, while the number of protons can be partitioned 
in many different ways, obtaining different results for the ratio defined in (l2()|), 


the relations from (21) on are only valid when the condition T^ = F^ = r°/2 


After finishing this work, we learned that BufalinojSSone of the coauthors of Ref. (231 has proposed 
evaluating partition energies as half the Q-values in the 150 — 166 MeV range. In this way, she 
obtains good agreement with £v™ UDA for A from 5 to 9, whereas for A = 13,15, and 16 there 
is a 2 a discrepancy. Note that the last range for the Q-values implies unrealistically small proton 
separation energies. 
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2 4 6 8 10 12 14 16 18 

Mass Number A 

Fig. 5. (Color online) Ratio R , given by as a function of A for the partition energies: (a) 
and (b) E£ en . 


is fulfilled, which does not occur for E^ x . There is, however, always an energy 
l?“ n for which this condition is fulfilled, and which are listed in the last column 
of Table [4j The values of new R, displayed in the lower panel of Fig. [5] are not 
very different from the previous values shown in the same figure. However, the 
value of the parameter a and, consequently, the ratio r^/Tp, is quite different now 
resembling those obtained by FINUDA, as can be seen from row C in Table [5] This 
agreement is somewhat surprising and we can not draw any conclusion from it. One 
has just learned that: (i) a relatively small modification of the partition energies 
(from E^ ax to El xBa ) can lead to significantly different results for rS^/Tp, and (ii) 
a relatively sizeable modification of the partition energies (from E™J’ OA to E*^ en ) 
can lead to quite similar results for T^/Tp. 

It would be valuable to find the physical meaning of the fitting parameter b, 
which is different in the three cases discussed above. In this regard, how to arrive at 
(22) from (21) is not a trivial issue. One possibility is to neglect the last term in the 
denominator of (211, arguing, as was done in Ref.that the FSIs tend to remove 
protons from the high energy part of the spectrum (N^ SJ < 0) while filling the 
low energy region (N p SI > 0), with the net result that N F gi = N p SI + N p SI = 0. 


Therefore 


R = a 


n fsj 


(i + r°/ro)NO’ 
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which, when compared with (22), yields 


N< s/ - MN“ ( 1 + J| 


bAN° p , 


( 25 ) 


since the factor 1 + T^/Tp < 1.5 is unsubstantial for a qualitative discussion. At first 
glance, the last equation appears reasonable because the effect of the FSIs should 
increase with A. But, since b = 0.009 ± 0.013, it turns out that Np SI fN p ~ 0.01A. 
Such a small amount of FSIs looks unrealistic. It is even more difficult to interpret 
physically the negative values of b that we obtain in cases B and C. Evidently, the 
fact that the IPSM is unable to reproduce the low-energy spectra in no way affects 
the previous discussion of r^/Fp. 

A different derivation of the 2N branching ratio has been done at FINUDA 
quite recently,^ based on the analysis of the (TT~,p,n) triple coincidence events, 


and a y 2 fit similar to (22). The new result, r§/F^ M = 021 ± 0.10, is consistent 


within the errors with the previously obtained valued as well as with our result. 
Only the sum of events from all hypernuclear species are exhibited in this work, 
without presenting data for individual proton spectra, which would allow us to do 
a reanalysis similar to that done above. 


6. Summary and Conclusions 

Proton kinetic energy spectra of ((He, ^Li, ® Be, ” B, ^ 2 C, \ 3 C, ^ 5 N and \ 6 0, mea¬ 
sured by FINUDA a few years agowere evaluated theoretically for the first time. 
We conclude that, in all the cases, the magnitudes of the spectra strongly depend on 
the parameterization that is used for the transition potential, while their shapes are 
very similar and independent of the transition mechanism. This statement, like all 
statements made in this work, do not depend at all on the inclusion or non-inclusion 
of the FSIs and 27V-NM channel in the nuclear model. 

It is explained in detail in Sec.[2]that our method for comparing the theory with 
experiment is radically different from the procedure followed by other researchers. 
In particular: 


The equations Q have never been used so far by any other group. These 
relations show that, to evaluate the decay rates, it is essential to know the 
number of produced hypernuclei N^, which is not available in the literature. 
We focus our attention on measured transition probabilities Tjv, and T n N, 
instead of comparing bare quantities F°, , r° p , etc , which are extracted from 

the experiments by making assumptions and approximations that are often 
questionable. 

The difference with other studies can be immediately seen by facing our Eq. 


(14) with 36 Eq. (7)]. 


But, what is really important is that the method proposed here imposes more 
constraint in comparing theory with data, allowing us to examine more clearly 
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the decay mechanism regardless of the importance of FSIs and the 2N-NM 
channel. 

Despite the lack of direct information about the N^-values, we have been able to 
estimate these observables for |He, ^Li, and ^ 2 C hypernuclei from the ratios R p 
presented in Refill (The physical meaning of this ratio is also clarified.) In this 
way, we obtain in Sec. [3] some very useful information on the transition potential. 
In fact, from Figs. [2] and [3] and Table [l] we conclude that: 

• The IPSM reproduces correctly the shapes of all proton kinetic energy spectra 
S p (E p ) at medium and high energies (E p > 40 MeV), including the Q- value, 
which is around 115 — 135 MeV. The latter could indicate that the residual 
nucleus A ~ 2 (Z — 1) is emitted mainly in the ground state. 

• This simple model also reproduces fairly well the magnitudes of the spectra at 
these energies when the soft tt + K potential is used to describe the NMWD. 
In no way do we claim that this is the ’’true” physics, but we believe that it 
might be worth pursuing this direction, especially considering that the model 
is able to explain satisfactorily the NM decay rates and T n of the s-slrell 
hypernucleiHE 

• The measured transition rates, r p , the same as the corresponding spectral den¬ 
sities Sp(Ep) for E p > 40 MeV, are significantly overestimated by the present 
theoretical calculations when the standard parametrization for the transition 
potential is used. In doing this comparison, one should keep in mind that, while 
the calculations refer to the 1V-NM decay mode only, the measured rates in¬ 
clude also the 2IV-NM decay channel and the effects of FSIs, and therefore the 
latter should always be larger. This happens only for the soft 7r + K exchange 
potential. 

• It is difficult to reconcile the FINUDA data with the theory based on the 
7T + 27 r/p + 27r/cr + io + K + pir/ai + air/ai exchange potential. 17 

We strongly believe that, in recent theoretical calculationswhich include 
both the 2V-NM decay channel, and the FSIs, they would have arrived at very 
similar conclusions if the comparison with experimental data have been made in 
the manner proposed here. 

Since the calculated spectra only involve the 17V-NM channel, the differences 
between them and the experimental spectra, both with respect to the magnitude 
and in relation to the energy distribution, can indicate which other degrees of free¬ 
dom are important (2V-NM channel, FSIs, etc ). Our plan of action is to add them 
to 17V-NM, when necessary, in order to successfully reproduce the experiments. 

With regard to the FINUDA method^ to determine the 2V-NM decay rate, 
based on the assumption that the 1IV-NM strength of the kinetic proton spectra is 
equally distributed in the low and high energies, we conclude that: 

• The proposed method is very sensitive to the energies that separate these two 
regions, and these energies can not be determined experimentally. 
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The separation done in Ref. 24 -'is not supported by any firm physical argument. 
It is necessary to resort to theoretical models to establish the partition energies; 
the IPSM is very suitable for this purpose. 

Both theoretical sets of partition energies (E^ h ax , E^ en ) differ from the FINUDA 
result (-E™ UDA ), not only in magnitudes but also with regards to the mass- 
number dependence: E^ x and £ , ‘ xen decrease with A, because the experimental 
single-particle energies increase; meanwhile, no reasonable explanation exists 
for the opposite behavior of E™ UDA . 

In spite of the above mentioned differences, all three sets of partition energies 
yield similar results for the parameter a and therefore for the ratio T^/Tp. This 
indicates that the behavior with A of the partition energies does not play a 
crucial role, and is consistent with a recent proposal to approximate them by 
a rather constant value of 1$ 80 MeV P® 

Physical interpretation of the FINUDA parameter b is a point at issue, not only 
for being very small, but also because it is negative in our analysis, as well as 
in the a new stud}® of the contribution of the 27V-NM channel employing the 
same method. Although suffering from large errors, its small value inevitably 
leads to the conclusion that the FSIs are very small, and that, therefore, the 
low-energy proton spectra dominantly comes from the 2IV-NM decay. It is very 
hard to understand this fact, and the only alternative possibility is that the 


basic approach (221 was incorrect. 


Final remarks: 


(1) As stated in the beginning, our purpose was not to reproduce the experimental 
data, but to discover out what the proton kinetic energy spectra can tell us 
about the weak hypernuclear interaction. To account for the low energy data 
of the kinetic energy spectra, it is imperative to consider the FSIs. For the time 
being, we are working on this issue by employing an improved version of the 
CRISP internuclear cascade model® used previously to describe C ITilLllJ 
Moreover, a complete theoretical description must also include a judicious es¬ 
timate of the effect of the 2./V-NM decay channel. So far, this has been done 
only in the context of FGMjHEH and it would be interesting to know what 
the SM can tell us about this process. In fact, for quite some time, we have 
been involved in the development of a corresponding theoretical formalism and 
numerical codes!® 

(2) To study the FSIs and the 2.ZV-NM decay, it is indispensable to understand first 
the liV-NM-decay dynamics. In our opinion, the SM could be a very useful tool 
to achieve this goal. For instance, the SM spectra are the main ingredients for 
establishing the initial conditions for the FSIs within the many-body multi¬ 
collision Monte Carlo cascade schemeEEH On the other hand, from Fig. |T| it 
is self-evident that the diagram a) plays the principal role within the diagram 
b) representing the 2A^-NM decay mode. 
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(3) New experimental developments will be very welcome, such as: i) Angular corre¬ 
lation of np and nn pairs to determine the and T nn rates, which have been 
measured so far only by KEK:® 4 ! in ^ 2 C, and ii) Triple (p, n,n), and ( p,p,n ) 
coincidence detections for direct measurement of F pn , and r pp , as suggested 
previously.^ The first steps in this direction seems to be given recently in 
RefP 

(4) To complete Figs. [ 2 ] and [ 3 J we need the Np/-values for ® Be, B, \ 3 C, \ 5 N 
and ^ 6 0 hypernuclei. Hopefully, these numbers will soon be available for public 
use. Needless to point out that, otherwise, the proton kinetic energy spectra 
measured by FINUDA^I are of little use to study the NMWD dynamics. They 
only can be exploited to discuss the decay kinematics through the analysis of 
the FSIs and the 2IV-NM decay mode. 
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